Fluctuation Study of the Specific Heat of Mg^^B2 
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The specific heat of poly crystalline Mg^^B2 has been measured with high resolution ac calorimetry 
from 5 to 45 K at constant magnetic fields. The excess specific heat above Tc is discussed in terms 
of Gaussian fiuctuations and suggests that Mg^^B2 is a bulk superconductor with Ginzburg-Landau 
coherence length = 26 A. The transition- width broadening in field is treated in terms of lowest- 
Landau-level (LLL) fiuctuations. That analysis requires that = 20 A. The underestimate of the 
coherence length in field, along with deviations from 3D LLL predictions, suggest that there is an 
infiuence from the anisotropy of Bc2 between the c-axis and the a-b plane. 
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Experimental observations of thermodynamic fluctua- 
tions in the speciflc heat have been limited in low-Tc su- 
perconductors because the long coherence lengths make 
the excess-Specific heat very small compared to the mean- 
field ternxl. By contrast, the high transition tempera- 
tures and small coherence lengths of cupratc^supercon- 
ductors lead to significant fluctuation effecjbffi. In the 
recently discovered superconductor Mg^^B2a, the coher- 
ence length and superconducting transition temperature 
lie between these extremes, suggesting that fluctuation 
effects will be observable and lead to further informa- 
tion on the superconducting coherence length. Indeed, 
the excess magnetoconductance of Mg^^B2 was reports 
recently and discussed in terms of fluctuation effects! 
Here we report the speciflc heat of Mg^^B2 from 5 K to 
45 K at several magnetic flelds. Using high resolution ac 
calorimetry, we could study the superconducting transi- 
tion region in detail. At zero-fleld, the excess specific 
heat is treated in terms of 3D Gaussian fluctuations and 
in flcld, the broadening and shift of the transition is ana- 
lyzed in terms of lowest-Landau-level (LLL) fluctuations. 

Polycrystalline Mg^^B2 was prepared at T = 950 C 
and p — Z GPa from a stoichiometric mixture of Mg 
and '^^B isotope using a high-pressure synthesis method. 
Since the sample was synthesized at high pressure, there 
has been no additional anaealing. Details of the synthesis 
can be found elsewhereO'Q. 

Measurements of the heat capacity were based on an 
ac-calorimetric techniqueQ. A long cylindrical sample was 
cut into a disk by a diamond saw and then was sanded to 
a thin rectangular shape whose dimensions are approxi- 
mately 1.1 X 1.5 X 0.1 mm'^; its mass is 375 fig. The front 
face of the prepared sample was coated with colloidal 
graphite suspension (DAG) thinned with isopropyl alco- 
hol to prevent a possible change of the optical absorption 
properties of the sample with temperature. The sample 
was weakly coupled to the heat bath through helium gas 
and suspending thermocouple wires. As a heating source, 
we used square- wave modulated laser pulse. The oscillat- 
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FIG. 1: Temperature dependence of AC at zero field, 
dash-dotted fine is a BCS fit with ACcxp/7n7c = 0.7. 
star represents a 3D Gaussian fiuctuation model above Tc of 
38.4 K. Inset: The temperature dependence of the specific 
heat at and 7 T from 5 to 45 K. 



ing heat input incurred a steady temperature offset (or 
dc offset) from the heat bath with an oscillating temper- 
ature superposed. The ac part was kept less than 1/10 
of the dc offset and was then converted to heat capac- 
ity by the relationship : C a I/Tac- The heat capacity 
obtained was converted to a speciflc heat by using a li|tp 
erature value above the superconducting temperatureB. 
The frequency of the periodic heating was chosen so that 
the ac temperature was inversely proportional to the fre- 
quency and, therefore, to the heat capacity; 23 Hz was 
used in this experiment. The ac and dc temperatures 
were measured by type E thermocouple, which were var- 
nished on the back face of the sample using a minute 
amount of GE7031 diluted with a solvent of methanol 
and toluene. The GE varnish typically amounts to less 
than 1 % of the sample mass. Since the field induced er- 
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FIG. 2: The excess specific heat, C/; = ACoxp(O), for T > 
Tc is plotted against reduced temperature t — T/Tc — 1 on a 
log-log scale. The solid line describes Gaussian fluctuations 
for a bulk superconductor with a G-L coherence length of 26 

A. 



ror of type E thermocouple is less than 1 % at 40 K in 8 
tesla, we will neglect the field dependence of the addenda 
contribution (DAG, GE-varnish, and type E thermocou- 
ple) and treat the field dependence as due only to the 
sample. 

The inset in Fig.l shows the temperature dependence 
of the specific heat at zero and 7 tesla from 5 to 45 K. 
The main graph is a plot of ACcxp vs temperature at 
zero field. Here ACcxp is the measured difference be- 
tween the mixed- and normal-state specific heats. A 7- 
Tesla data set was used as a reference state above 20 K 
because it shows no observable transition in that range. 
The subtraction was executed without any smoothing of 
the 7-T data. The dash-dotted line k a BCS fit with the 
ratio of ACcxp/"/nTc being variables. The normal elec- 
tronic coefficient 7„ was set to be 2.6 mJ/mol K from the 
literatureH and Tc of 38.4 K was determined from scaling 
discussed below. The best fit showed that the ratio is 
0.7, which is much smaller than the weak coupling BCS 
value of 1.43. Since the ratio is generally larger for strong 
coupling superconductors, the small value does not tell 
us anything about its coupling strength. Recently, there 
has been a plethora of experimental and theoretical ev- 
idence which supports two-gap features in MgB2, which 
can ex]jJaiji pth^ Jipn-BCS jump magnitude with some 
success Jl3'U't3ll§llJ However, we cannot rule put— such 
other scenarios as an anisotropic gap structureJl3'E£l For 
a system in which fluctuation effects are pronounced, 
the experimentally determined transition temperature is 
lower than the mean- field critical temperature (T^^) be- 
cause fluctuations drive the system into the normal state 
even below T™-^ . It is unlikely, however, that this can 
explain the large deviation from the BCS value. 

Above the transition temperature, there is an ept. 
cess specific heat tail apparent in Fig. 1. ThoulessEII 
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FIG. 3: AC'/T is plotted against temperature in B=0, 0.2, 
0.4, 0.6, 0.8, 1.0, 2.0, and 3.0 tesla. Data were taken with 
increasing temperature after field cooling. The solid line rep- 
resents a BCS fit. 

and subsequently Aslamazov and Larkin0 showed that 
Gaussian fluctuations arise above Tc and predicted that 
Cfi = C+<-{2-'i/2) ^iti^ (j+ ^ {kB/STT)(GLiO)-^, where 
t = T/Tc—1, d is the dimensionality, and 'Cgl(O), the T=0 
K Ginzburg-Landau coherence length. Figure 2 shows 
the temperature dependence of the excess specific heat 
on a log-log scale. The data follow a power law with 
an exponent of -0.5 and C+ — 0.66 mJ/mol K. The ex- 
ponent indicates that Mg^-'^B2 is a 3D superconductor 
and the substitution of C^ into the above formula gives 
fGL(O) = 26±1 A. 

When a magnetic field is applied, the specific heat 
broadens. Figure 3 shows the temperature dependence of 
AC/T at several magnetic fields. The ratio of the tran- 
sition temperature shift to the transition width broad- 
ening in field is unique in that it is not as large as in 
low-Tc superconductors nor as small as in high-Tc ma- 
terials. Its intermediate behavior is related to the fact 
that the coherence length and the superconducting tran- 
sition temperature of Mg^^B2 are intermediate betweeu. 
low-Tc and high-Tc superconductors. Lee and Shenoyu2l 
studied fluctuation phenomena in the presence of a mag- 
netic fleld, arguing that bulk superconductors exhibit a 
field-induced effective change to one-dimensional behav- 
ior in the vicinity of the transition temperature Tc{B). 
In a uniform magnetic field, the fluctuating Cooper pairs 
move in quantized Landau orbits and, close to upper crit- 
ical field (Bc2), the lowest Landau level dominates the 
contribution to the excess specific heat. So, a bulk su- 
perconductor behaves like an arraK-.of one-dimensional 
rods parallel to the field. ThoulesaHJ extended the idea 
above and below Tc and suggested a scaling parameter 
for the fluctuation specific heat that is valid throughout 
the transition region : 




3 



< 



0.6 ■ 




B{T) 




T 





38.4 


0.008 


0.2 


37.2 


0.009 


0.4 


36.5 


0.012 


0.6 


35.5 


0.0165 


0.8 


34.8 


0.02 


1.0 


34.3 


0.025 


2.0 


32.5 


0.042 



t/x 



FIG. 4: ACcxp/ACm/t vs t/r is shown, where ACcxp ~ 
C{B) — C{7T) and ACmft are defined in the inset. Here t is 
reduced temperature and r is a field dependent dimensionless 
broadening parameter. Sohd squares are from a Hartree-like 
approximation. Other field data were scaled so that they were 
collapsed onto the Hartree result. The obtained parameters 
(Tc, r) are tabulated inside the figure at each magnetic field. 
The inset sketches the definition of ACmft{T, B)E3 



where t is the reduced temperature and t is a field depen- 
dent dimensionless parameter that describes the super- 
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conducting transition width. The functional form g{y) 
model dependent. When a Hartree-like approximatior 
is used to examine the fluctuation efi^ects of the quartic 
term in the free energy functional, it results in a simple 
form : 



where y — t/r^B). 

In Fig. 4, the ratio of ACcxp/ ^Cmft 
region was plotted as t/T{B), where AC, 
C(7T) and AC™/* wai 
by Farrant and Gougl 



(2) 
(3) 

in the transition 
cxp(B) =C(B)- 



.etermined as in the classic work 
I by fitting the low temperature 



side of ACcxp (S) as in Figs 1 and 3 and extrapolating 
linearly above Tc{B). A sketch is shown in the inset 
of Fig. 4. The scaling parameters t{B) and Tc were 
chosen to make the data collapse onto the Hartree-like 
approximation (solid-squares). The values of t{B) and 
Tc are listed in Fig. 4. The temperature dependence of 
the upper critical field Tc(B) is plotted in Fig. 5, and 
shows positive cupjature close to Tc{B = 0). A simple 
empirical formula^, Bc2(T) = Bc2(0)[l - {T/Tc)^][l - 
a{T/Tc)^], was used to describe the curvature, in which a 
is a fitting parameter ihat is and 0.3 for two-fluid model 
and for WHH modeO respectively. The best fit, solid 
line in Fig. 5, was produced with i?c2(0) = 15.4 tesla 
and a = 0.8. Positive curvature near Tc{B = 0) was also 
observed in non-magnetic rare-earth nickel borocarbides 
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FIG. 5: The temperature dependence of the upper critical 
field Bc2(r) is shown. The solid line is a simple empirical 
formula, Bc2(0)[l - {T/nf][l ~ a{T /T^f], with a = 0.8 and 
Bc2(0) = 15.4 T. The dotted line is a linear fit with Bc2(0) = 
16 T and dBc2/dT = -0.44 T/K. 



RNi2B2C{R = Lu,Y) and could be explained by the 
dispersion ot|thc Fermi velocity using an effective two- 
band modelcJ. 

The broadening parameter t(B) consists of a field de- 
pendent part (tb) and a field independent part (Ti„). 
We postulate that they are independent of each other 
and add in quadrature, such that the total broadening 
is = + T^. The field independent part was ob- 
tained in zero field and accounts for sample inhomogene- 
ity and zero-field fluctuation effects while the field depen- 
dent part is due solely to field-induced fluctuations. The 
field, dependence of the broadening parameter is given 
hyM: 
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where a = 2 — (d — 2)/2 {a = 3/2 for a bulk super- 
conductor) and Bg ~ ~Tc{dBc2/dT)x^ in the mean-field 
scheme. The exponent d — 2 indicates a dimensional 
crossover from d-dimension to d-2 behavior. The shift 
field Bg is a characteristic field that sets the scale of 
the shift of the transition temperature while Bw sets the 
scale of the width broadening of the transition region. 
In a standard superconductor, the ratio By/ 1 Bg is very 
large 10'*), and the transition is shifted much more 
rapidly in field than it is broadened. In high temper- 
ature superconductors, such as YBCO, the broadening 
is as large as the shift {Bg ~ Bw), which is an indi- 
cation that a mean-field approach based on a perturba- 
tion expansion might not be proper and that fluctuations 
should be treated in the context of critical phenomena. 
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FIG. 6: Field-broadened dimensionless parameter tb is plot- 
ted against B/AC on a log- log scale. The parameter was ob- 
tained by using tb ~ (t^ — rj^„)^''^ and AC is a mean-field dis- 
continuity at Tc. The solid line represents a lowest-Landau- 
level (LLL) fluctuation: tb = /3(-B/AC)^/" with a = 3/2 and 
/3 = 2.7 [A/m-Kf/^. 

In Mg^'^B2, the ratio is in the order of 10^, which is in be- 
tween those two extremes. This feature seems consistent 
with other properties that show aspects of both conven- 
tional and high Tc superconductors. In order to study 
the anomalous broadening in field, we plot tb vs (B/ AC) 
on a log-log scale in Fig. 6. The slope represents the 
exponent (l/a) while the coefficient of the slope (3 is re- 
lated to the Ginzburg-Landau coherence length through 
(3 = (fcs/STT^giJs)^/". The lowest-Landau-level approxi- 
mation is shown as a solid line having a slope of 2/3 and 
coefficient 2.7 (A/m K)^/^. From the coefficient of the 
fit, the Ginzburg-Landau coherence length is estimated 
to be 20 A. In the above analysis, the shift field Bs of 16 
tesla was obtained by fitting the linear region of Bc2 (see 
Fig. 5). ^ 

Before making any further conclusions, we discuss 
some of the assumptions that we made in the analysis. 
In the zero-field analysis, sample inhomogeneity has been 
neglected. Boron 11 isotope Mg^^B2 has a Tc of 39 K 
while the excess specific heat extends well above 40 K. 
We might expect inhomogeneity effects to complicate the 
fluctuation analysis below the transition temperature but 
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above Tc, where our analysis is concerned, the effect will 
be negligible. However, when nonzero field is applied, 
sample inhomogeneity must be considered because the 
analysis is of the field dependent behavior of the tran- 
sition region. Sample inhomogeneity could produce an 
additional broadening through the Ginzburg-Landau pa- 
rameter K , and hence Hc2. Information on Hc2 slopes 
at different parts of the sample with different Tc's would 
be needed to account for the additional broadening cor- 
rectly. To be more precise, the Hc2 slope in the Tc = 39 K 
part of the sample and that in the, e.g., Tc = 38 K part 
would be needed. In our in-field analysis, we assumed 
that the slopes of Hc2 at different parts of the sample are 
same or if they are different, the difference is small, which 
leads to field independent inhomogeneity effect. It is nec- 
essary to study high quality single crystals with different 
Tc's to better understand sample inhomogeneity effects 
on transition-width broadening. 

In summary, the zero-field specific heat was discussed 
in the context of BCS theory plus 3D Gaussian fluctua- 
tions. The analysis indicates that Mg^^B2 is a bulk su- 
perconductor and its coherence length is about 26 A. In- 
field specific heat was treated in terms of lowest-Landau- 
level fluctuations. That analysis requires that <^o — 20 
A. The in-field analysis could be complicated due to the 
effect of anisotropy in a polycrystalline sample. The 
anisotropy of i?c2 between ab-plane and c-axis directions 
can lead to a field dependent broadening due to the Tc 
distribution arising from the randomly oriented grains of 
the present sample. This in turn leads to an underesti- 
mation of the G-L coherence length. In order to under- 
stand the influence of anisotropy, we can assume that the 
transition broadening arises solely from _Bc2 anisotropy 
and calculate the required anisotropy in our experimental 
temperature range. The reportecLamsotropy value of ^3 
from single crystal measurementsc3^E3'E3 is much smaller 
than the ratio 6 required to explain the broadening. From 
this consideration, we conclude that the anisotropy alone 
cannot explain the whole broadening and therefore that 
fluctuation effects should be considered in explaining the 
anomalous broadening. 
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